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Abstract 

We apply the thermodynamical model of the cosmological event hori- 
zon of the spatially flat FLRW metrics to the study of the recent acceler- 
ated expansion phase and to the coincidence problem. This model, called 
"ehT model" hereafter, led to a dark energy (DE) density A varying as 
r~ 2 , where r is the proper radius of the event horizon. Recently, an- 
other model motivated by the holographic principle gave an independent 
justification of the same relation between A and r. We probe the theoret- 
ical results of the ehT model with respect to the Snla observations and we 
compare it to the model deduced from the holographic principle, which we 
call "LHG model" in the following. Our results are in excellent agreement 
with the observations for Ho = 64fcms _1 Mpc~ 1 , and Q.^ — 0.63lo Jii 
which leads to go = —0.445 and zt — 0.965. 
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1 Introduction. 

Since the discovery of the presently accelereted expansion of the universe from 
supernovae observations £P[2], evidences for such an accelerated phase are in- 
creasing. The simplest theoretical candidate to explain this acceleration is a 
cosmological constant A. Anything producing sufficient negative pressure - for 
instance a scalar field [3] or a bulk viscosity 0] - could also be valid. 

Before the discovery of this acceleration, phenomenological ansatze with a 
variable A(i) were tentatively proposed as solutions of the cosmological "con- 
stant" problem. For instance, laws such as A(t) ~ t~ 2 A(t) ~ a~ 2 (t) HUH, 
where a(t) is the scale factor of the FLRW space-time, A(t) ~ H 2 (t) 8 or 
A(t) ~ (3H 2 {i) + (1 - / 3)iJ 3 (t)i?7 1 9 , where H(t) is the Hubble parameter and 
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Hi is its exponential inflation value, were suggested. Other proposals |1 | |1 1 | 
can also be quoted. 

From a different point of view, the generalization J2| EH of the black hole 
and of the de Sitter event-horizon Thermodynamics |14U15| to the FLRW space- 
time has led to the relation A(t) ~ r~ 2 (t) ^H] where r denotes the event-horizon 
in the FLRW model of the universe. 

Recently, this last form for A(t), or, equivalently, for the dark energy den- 
sity p\{t) through p\(t) = x _1 A(£) (with \ = 87rGc~ 4 ), has received further 
confirmations based on the holographic principle ^| EH • 

Let us remind of the approach which can been followed to produce a model 
with a time-dependent cosmological constant. We start with a type-like perfect 
fluid energy-momentum tensor 

T al3 = ptot^u 13 - P tot A al3 , A af3 = g aP - u a u^ (1) 

where u a is the 4-velocity common to all the components of the energy 
density p tot . We consider two components such as p to t — p + PA and P to t — 
P + P\. The component {p, P) is the matter with the barotropic state equation 
P = (7 — l)p where 7 is a constant (for instance, 7 = 1 for dust). The second 
component is the Dark Energy (DE) with p^, the vacuum energy density, and 
Pa the (negative) pressure satisfying the state equation 

Pa = -pa- (2) 
Relation (2) leads to the two following alternatives: 

i) Each component is conserved separately and, of course, A has to be con- 
stant. 

ii) Both of the components are conserved together, A = A(i) is then possible. 

The event-horizon Thermodynamics (ehT) model is derived on the basis of 
point ii) by assuming an interaction between the matter and the Dark Energy 
(DE hereafter). Let us remark that we write "matter" for any sort of matter 
except DE. Today, the matter is the dust, the largest part of which is the 
Dark Matter (DM). For sake of simplicity, we use DM to denote the dust, 
encompassing the baryonic matter. In the same vein, other models assuming 
an interaction between the DE and DM components of the cosmic fluid were 
studied, e.g. [T§|. 

A model such that A ~ r~ 2 for the DE density can be used in different ways 
and different contexts. For instance, in a precedent paper |16| in order to address 
the problem of the exit of inflation in the early universe, we imposed as second 
component a perfect fluid of strings (7 = 2/3). The model led then to A = 
3^, which was independently considered as an ansatz derived by dimensional 
considerations by some authors [2U| |2J An equivalence can be found 

between the previous relation A ~ - and the forms A ~ a~ 2 and A ~ p under 
specific conditions |23| . 
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In the present paper, in order to settle some issues on the coincidence and the 
recent decceleration-acceleration transition problems, we assume for the second 
component a cold dark matter (P = 0). In section 2 we review some basic 
equations and relations common to the ehT and LHG models. The ehT model 
is developed in section 3, particularly for the z < 2 epoch. In section 4, in order 
to probe the DE assumption in this range of z, we discuss how our model fits in 
with the type la supernovae observations [21] ■ We deduce then the most likely 
values for the Hq and f2 A parameters, as well as the deceleration parameter qo 
and the decceleration-acceleration transition redshift zt- Finally, sections 5 and 
6 contain comments and a brief comparative discussion concerning the results 
obtained by the two models. 



2 Model for A and Field equations. 

In order to set the notations, we introduce some basic equations of the two 
component models. The spatially flat FLRW space-time has the metric 

ds 2 = c 2 dt 2 - a 2 (t)[dR 2 + R 2 (d9 2 + sin 2 6d(j> 2 )} (3) 

where the scale factor a(t) is a monotonic increasing function of the cosmic time 
t. 

We assume an universe filled by two interacting type-like perfect fluids, 
namely dust (ordinary and dark matter) and Dark Energy (DE) . The dust and 
DE energy densities are p and pa = X _1 A respectively, and their correspond- 
ing pressures are P and Pa. The two state equations are P = (7 — l)p with 
7 = const, < 7 < 2 , and Pa = Lup\, where u) can be variable. 

We recall the field equations for the spatially flat case 

3P 2 = X c 2 {p + PK) (4) 



2- + P 2 = -Xc 2 (P + P A) , (5) 

a 

where H = - is the Hubble parameter, c the velocity of the light and the dot 
stands for the time derivative. 

Combining these two equations leads to 

(ff'- 1 ) = |( 7 +(l + ^-7)^A) , (6) 

where the dimensionless density parameter £7a = Ac 2 /3H 2 has been introduced. 
The equation (6) is always valid provided the DE is a perfect fluid. 

We consider now A as a vacuum energy density associated to the FLRW 
event-horizon such as 

a 3q;2 
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where r is the proper radius of the event-horizon, and a is a dimensionless 
constant parameter. This form of A was previously obtained by and [T7| 
when a = 1, and by |18| when a^l. 

Using the quantity ft a, relation 0) becomes 

The proper radius of the flat FLRW event-horizon is 

r ( t } = <t )f^ y (9) 
The derivative of J5J) with respect to time gives 

T C 

H-- = -. (10) 
r r 

For convenience, we introduce the variable x = lna(i) such as x = today. 
Relation (|l(Jfl becomes then 

1-^ = ^ , = (ID 

where the prime means the derivative with respect to x. 
In the same manner, we can rewrite relation (JBJ 

7^7 = 1(7 +(1+^-7)00- (12) 

Finally, by combining equations (jllfl and l|12fl with the derivative of equation 
JSJl, one obtains 

A = f> A {3[7+(l+^- 7 )f! A ]-2[l-^]} . (13) 

a 

Let us emphasize that this equation is valid for any values of 7 (constant) and 
u> (constant or variable) , independently of the fact that the two components p 
and pa are interacting or not. 

It is useful to derive from the field equations J3J and (JSJ the deceleration 
parameter q 

«=^=i[(37-2) + 3(a; + l-7)fiA]. (14) 

which is valid in the two models. 

In the following, we assume that the "matter" component p is dust (7 = 1), 
so that l|13|) and i|14|l become 

n' A = n A (i + 2^—^ + 3ojn A ) (is) 
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g=-(l + 3wfi A ). (16) 

The relations il3*|)- (|16[l are valid in the two models under consideration, which 
we denote K(t)CDM models hereafter. 

From now on, the assumptions of the ehTmodel will be different from the 
LHG model's ones. 



3 Model with interacting components. 

We assume that the DE component satisfies thermodynamical state equations, 
i.e. relations between its thermodynamical variables which are valid in any 
space-time. Therefore, any thermodynamical state equation valid in the de 
Sitter's space-time ^HIEH] - f° r instance, Fa = — PA and pa — 127r 2 T A 2 (T a 
the temperature ) - remains valid in the FLRW space-time. Thus, if the DE 
is an actual cosmological component, its thermodynamical state equations will 
stay the same, independently on the choice of the space-time as well as for 
any other component. This suggests to retain the relation J7J) which is valid 
in the de Sitter's space-time when a = 1. In section 5, some consequences of 
the presence of the parameter a in the ehT and LHG models are discussed. 
Using the holographic principle can lead also to choose the relation (JJJ [T7] . 
[T%] . These references assume a variable state equation (oj — oj{x)) for the DE, 
and independent energy conservation laws for the matter and DE components. 
Conversely, the present model assumes oj = — 1 (vacuum), and that the energy 
conservation is only valid for the two components considered together. 
Equation (|15JI can be rewritten 

n' A = 3f> A (/3 2 - V^aXA + V^) (17) 
where the constants /3iand P2 are given by 

/?! = -!-(Vl + 3a 2 -l) , /3 2 = -L(v/l + 3a 2 + l), /?i,/3 2 >0. (18) 
6a 6a 

By setting a — I, Equation (17) becomes 

n A = n A (i- v / ^a")(3 v / ^a" + i), (is) 

which differs from Equation JSJ in . Nevertheless a straightforward calcula- 
tion (using ifT^l) . l(T5|l and the derivative of the definition of Ha) gives 

A' = 2A( V / Oa - 1) , (20) 

which is common to the two models. As A' is always negative, A is decreasing 
with time. Observational evidences provide a very small present value for pa 
(fine-tuning problem) and of the same order as p (coincidence problem). 
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Introducing the function y(x) = V^Ta, Relation l|17l) becomes 

2tf , = 3v0%-y)(/9i+ 1 /). (21) 
Its solution is (in the only case considered here where y < $2) 

Km= s y a (22) 

_ y) Vi+iSr ^ + y) ft 

i^i is a constant of integration which can be related to the initial condition 

yo = v^a- 

We derive now the expression of r = r(y). Using Equations (|llf> and l|21|l 
yields 

d (In r) = dx- - , a 2 % a r . (23) 



After integration, one obtains 



K 2 r = a(^-^)V^ (24) 
or equivalently 

Kr= a _, 2 y 2 a+fii , K^K X K 2 . (25) 

(#2 - y)f 2 Vi+3° 2 (fa + y)ftVl+3 Q 2 

2<Ta is a second constant of integration which depends on yo and ro — ac(yo^o) ■ 
The expressions of K\ and #2 depend explicitly on the two priors and Hq. 
The current values of Vi\ and i/o are = 0.7 and Hq = 72km.s^ 1 .Mpc 
|26| . With these two numerical values, it is interesting to deal with the case 
where a = 1 for which 0i = -g and (3 2 = 1. One obtains 

#ia = 5- , #1 = P-. 3- = 1.3686 (26) 

{i-y)H\+y)i (i-w>)*(£ + w>)* 



K 2 r = ad — - )5 , or Kr = (-r^— ) 2 , r = — — = 4980.12Mpc (27) 

g + y ^ + y H oyo 

K 2 = — ( ~ V ° )2 = 7.50265 x 10~ 5 Mpc' 1 , X = 1.02681 x lO^Mpc -1 . 

However the previous values of Hq and 17" are model-dependent. They were 
obtained in the framework of the ACDM model. We shall see that starting 
with the same observational Snla data, the best fit to the K(t)CDM models 
give appreciably different central values of H and f2° . 
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4 Snla constraints on the ehT model 



In order to compare these theoretical results with the observations of the Snla 
magnitudes, the luminosity distance dh has to be expressed with respect to the 
redshift z = a^ 1 — 1. In the ehT model, it yields 

d L = (l + z)[(l + z)r r„] - ^-^[(1 + z)~ - 1], (28) 

where the expression of r depends on 2. As before, we only consider the case 
a = 1. Both Equations (1221) et H25|) give a parametric representation (via the 
"parameter" y) of r as function of z. Indeed, (22) yields immediately z = z(y) 
(with a = (1 + z)- 1 ). 

The set of the theoretical curves "distance moduli" fi versus the redshift z, 

H = m-M = 25 + 5 log 10 (d L ) , with d L in Mpc, (29) 

predicted by the model parametrized by the two cosmological parameters yo = 
\/^a et ^Oj can be plotted. For the two parameters f2 A and H free, the best 
fit to the magnitude observational data of the 157 Snla "Gold sample" PU 
can be determined by minimizing the function X 2 = E C^"^' ) 2 , where 
fJ-i(zi) denotes the values of the magnitude for the observational data, Oi the 
corresponding error and the summation is taken over any of the 157 data of 
the sample. The corresponding values of and Hq are derived by numerical 
computation. More precisely, Equation l|21|) is integrated by the method of 
Runge-Kutta of order 4, and the expression of z{y) is deduced by use of 
With the help of Equations Qffi and J2SJ|, the values of fj,(z) for z ranging 
from to 100 are then obtained. After a simple numerical evaluation of \ 2 for 
fi A ranging from to 1 and Hq from 50 to 100, the best fit corresponding to 
X 2 = 178, 7 is obtained for 

H a = 64+lkm.s- 1 .Mpc' 1 ,Q° A = 0.63±g;J x , (30) 

The function ^i(z) is plotted in figure 1 for z ranging from to 2. 

The likelihood function C(fl\) (see figure 2) is derived by marginalization of 
Hq and furnishes the same value of the parameter Vvi . 

Finally, the decceleration parameter q can be expressed as a function of y in 
the ehT model ( for a = 1) from equation (16) ( with lo = —1) 

q = \{i-Zy 2 ). (31) 

In figure 3 the curve q(z) of the ehT model is plotted. Today the decceler- 
ation is go = —0.445 , and the decceleration-acceleration transition occured at 
z T ~ 0.965. 
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5 The event horizon and the parameter a. 

We examine here the influence of the parameter a on the limits of the proper 
radius r of the event horizon (eh) in the two models. First, let us consider the 
LHG model. 

By comparison with the relations (|22[) and (|25|) of the ehT model, the LHG 
model would lead to the relations ( a is given by (9) of ^H] an d r i n °t explicitly 
given, can be deduced from their eqs. (6) and (9)): 

y 2 (i + y)^ yo(i + yo)^ 

Y a = — — , a ± 2 , Y - 



(1 - y) 2+» (a + 2y)^ 1 (1 - y ) {a + 2y a ) 1 ^ 1 

(32) 

a y 2 (i + y)^{i~y)^2_ (33 ) 



Y*H ^T^W K (a + 2y)^ 

For a = 2, the LHG model requires to start again the calculation from the 
differential equation i|15|) which becomes: 

2y' = y(l~y)(l + y) 2 . (34) 

Its integration yields 

(1-2/0)3(1 + 2/0)3 y 2 8 1 1 

Vl (l-y)Hl + y)i XPl 34+y 1 + yo"' 

Then, 



2c(l - y Q ) 2 (l + y ) V *M—y ~ T+^- 
iJoyr 3 ^?/;? (l-y)i(l + y)i 



(36) 
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Figure 2: The likelihood function £(0^) of the ehTmodel 



We can see from (|33[) or (|36|1 that a tends to infinity when y tends to 1, for any 
values of a (positive, see (8)). But the behaviour of r differs because it depends 
on the parameter a, as it can be seen from (32) and (35). Three cases can be 
distinguished for the behaviour of r in the limit y —> 1: 



if a < 1 



(37) 



oo if a > 1 



(38) 



9 tf yT^Y 5 



2 (1 + 2y ) 2 
9 (1 + 2/o)yo 



c 

Hi 



if a = 1. (39) 



The first two cases ( i.e. r — > and r — ► 00 ) disagree with the holographic point 
of view, because they would prevent any cut-off (IR and UV respectively). In 
particular, the case a < 1 seems to be proscribed because it could not prevent 
the singularity formation and would correspond to the absence of black hole 
formation . 

The third case only (a = 1) corresponds to a de Sitter asymptotic limit. In 
Equation l|39|) , the index i of H means exponential "inflation" . Note that the 

= 1.06813 if we take y Q = VoTf. 



limit ^ depends only on y , and its value is : 



As r = 



Hoyo 



4980.12Afpc, r, ; is equal to 5319.42Mpc. The expression of r 



is formally the same in the two models and depends only on the choice of the 
observationnal priors Hq et yo- However, each model leading to slightly different 
adjustments of these parameters gives slightly different values of r$ and then. 
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Figure 3: The decceleration parameter q(z) of the ehTmodel 



In the case of the ehT model, for any arbitrary a, the same phenomenon 
appears and the value a = 2 does not necessitate a special study. In the limit 
y — ► 1, Equations l|26|) and (|27[) give 



a — ► co and r — > csi = ^(^) 2 = 5478.13Afpc if a = 1 (equivalently, /?2 = 1) 

When a > 1, /?2 < 1, then y —> P2 before reaching 1, and a —* 00 , while r — > 00 
for this asymptotical limit /?2 of j/. From the today observational evaluations, 



(3 2 has to be > V0.63 = 0.79, and so a < g^o^i = L78 - In tnc future, a 
range from 1 to 1.78 will become more and more narrow, tending to 1, as long 
as the equation (17) of the model, indicating a growth of J7a, remains valid. 

Thus, the case a — 1 appears to us as the most attractive. The corresponding 
de Sitter's limit is r, = 5478.13Mpc. It is a little greater than the limit of the 
LHG model (5319.42Mpc), which means a little weaker exponential inflation. 

6 Conclusion. 

We have seen that the form A ~ r~ 2 , clearly supported by the holographic 
principle, leads, in our study, to two somewhat different models, owing to the 
chosen energy conservation equation. In the ehT model, a — 1 and the best fit 
(xt = 1.14) to the Snla's data from the "gold" sample [21] gives us Ho = 64 
km.Mpc^ 1 .s _1 and 57° = 0.63. If a ^ 1 (as in the LHG model) it is worth 
observing that the a < 1 values are not very attractive because they lead to the 
singularity r — > when Q\ — ► 1. 

For the decceleration-acceleration transition epoch we find a redshift zt — 
0.96, a value slightly higher than the ones recently published (0.28 ^ zt ^ 



a — > 00 and r — > if a < 1 (equivalently, /?2 > 1) 



1 ,3 
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0.72)^21 [21] an d very sensitive to the 17° value. Comparing the values of 
the cosmological parameters in various models requires to discuss not only the 
choice of the parameter a but also the forms or relations taken for q(z) (for 
instance, q(z) = qo + qiz valid when z <§C 1), for w(z), or for cLl{z). Besides, in a 
given model, one has to take into account the energy conservation laws for DM 
and DE. In most cases, the authors assume an energy conservation law for each 
component separately. Here we have considered the more general situation of a 
global conservation of the whole energy and, necessarily, an interaction between 
DM and DE. Such an interaction could induce higher values for the transition 
rcdshift Zr, as noted by Amendola et al. for models with coupling |^|[2S|. Fu- 
ture observations in the high redshift range could allow to discriminate between 
theories with coupled components and theories with distinct conservation laws. 
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